JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOL10.1029/, 



The role of oxygen ions in the formation of a bifurcated current 
sheet in the magnetotail 

S. Dalena 1 , A. Greco 1 ' 2 , G. Zimbardo 1 ' 2 , P. Veltri 1 ' 2 

1 Dipartimento di Fisica, Univcrsita dclla Calabria, Rcndc (CS), ITALY 

2 Consorzio Nazionalc Intcrunivcrsitario per le Scienzc Fisichc dclla Materia (CNISM), Unita' di Cosenza, Pontc 
P. Bucci, Cubo 31 C, Rcndc (CS), ITALY 

Abstract. Cluster observations in the near-Earth magnetotail have shown that some- 
times the current sheet is bifurcated, i.e. it is divided in two layers. The influence of mag- 
netic turbulence on ion motion in this region is investigated by numerical simulation, 
taking into account the presence of both protons and oxygen ions. The magnetotail cur- 
rent sheet is modeled as a magnetic field reversal with a normal magnetic field compo- 
nent B n , plus a three-dimensional spectrum of magnetic fluctuations <5B, which repre- 
sents the observed magnetic turbulence. The dawn-dusk electric field E y is also included. 
A test particle simulation is performed using different values of (5B, F, y and injecting two 
different species of particles, + ions and protons. + ions can support the formation 
of a double current layer both in the absence and for large values of magnetic fluctu- 
ations (SB/B = 0.0 and 5B/B > 0.4, where B is the constant magnetic field in the 
magnetospheric lobes). 



1. Introduction 

The magnetotail current sheet, which separates the north- 
ern from the southern lobe, is one of the key regions of mag- 
netospheric physics. In the simplest 1-D approximation, it 
may be described by the Harris solution, an equilibrium so- 
lution of the Vlasov equation, in wich the current density 
maximum is near z = 0, where the magnetic field is equal to 
zero. Although spacecraft observations confirmed the Har- 
ris model as a zero-order description of the current sheet, 
later studies and observations by spacecraft revealed a more 
complex structure. Previously, current density profiles with 
a double peak were observed in the near-Earth tail by ISEE 1 
and 2 [Sergeev et al, 1993] and were deduced by a statistical 
analysis of Geotail data in the distant tail [Hoshino et al., 
1996]. More recently, bifurcated current sheets have been 
reported with Cluster data at 20 Re downtail [Nakamura 
et al, 2002; Runov et al., 2003b]. In some of the analyzed 
events, the formation of the current double layer was asso- 
ciated with magnetic reconnection, as shown by Runov et 
al. [2003a] and Asano et al. [2004]. Observations of 29 Au- 
gust 2001 and 26 September 2001, analyzed respectively by 
Runov et al. [2004] and Sergeev et al. [2003], reported on 
a bifurcated current sheet during apparent flapping motion 
of the plasma sheet, that lasted for 15 min. The authors 
concluded that the bifurcation was not associated with re- 
connection process but was the result of an aging process of 
a thin current sheet, in which non-adiabatic motion of ions 
results in a weaker current in the central plasma sheet, as 
demonstrated theoretically by Zelenyi et al. [2002]. How- 
ever, the physical origin of bifurcated current sheets remain 
controversial. 

Contemporary to the observations, many theoretical 
models have been proposed. Zelenyi et al. [2002, 2003] con- 
sidered the structure of a thin current sheet in presence of 
non-adiabatic particles, wich reduce the current at the cen- 
ter of the current sheet and the bifurcation has scale size 
of an ion gyroradius (hundreds of kilometers). Sitnov et al. 
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[2003] obtained bifurcated current sheets in the case of ion 
temperature anisotropy T± > Ty (where Tj_(7j|) is the ion 
temperature perpendicular (parallel) to the magnetic field), 
together with relatively flat current density profiles. Ricci et 
al. [2004] developed a kinetic simulation of a Harris current 
sheet in which the lower-hybrid drift instability caused the 
current sheet to bifurcate. Delcourt et al. [2004] examined 
the behavior of charged particles with adiabaticity parame- 
ter k—1 in a simple current sheet model and found that the 
non linear dynamics lead to bifurcated current sheets on the 
scale of ion gyroradius. 

In most previous studies, the presence of magnetic tur- 
bulence in the near and distant magnetotail was not consid- 
ered. Although the magnetic turbulence is stronger during 
active geomagnetic periods, it is non-negligible even dur- 
ing the quiet times. Usually, the turbulence is stronger in 
the center of plasma sheet [Bauer et al., 1995; Hoshino et 
al, 1996], while a more ordered magnetic configuration is 
found in the vicinity of the magnetospheric lobes. Veltri et 
al. [1998] investigated the effect of turbulent magnetic fields 
on the proton dynamics for the distant magnetotail, where 
the normal component of the unperturbed magnetic field, 
Bo z = B n , is statistically negligible (in the GSM coordi- 
nate system). By a test particle simulation, they found that 
the magnetic fluctuations play the role of effective scatter- 
ing mechanism and the current sheet splits in two layers for 
perturbation levels SB /Bo > 0.2. Subsequently, Greco et 
al. [2002] and Zimbardo et al. [2003], analyzing the near- 
Earth magnetotail, where the average B n is not zero and 
northward oriented, found that the normal magnetic field 
component B n and magnetic fluctuations SB /Bo have op- 
posite effects on the current structure and on the proton 
heating. Indeed, a large value of B n inhibits the i/-motion in 
the quasi- neutral sheet, while magnetic fluctuations favour 
the y motion at some distance from the center; the strong 
magnetic turbulence in the center of the current sheet slows 
down the proton motion and causes the double humped pro- 
files of the current and ion velocity. More recently, Greco et 
al. [2007] have studied the equilibrium of the Earth magne- 
totail, injecting protons and electrons. A stationary three- 
dimensional kinetic-fluid code with protons represented by 
particles and electrons by a mass less fluid was developed. 
For a specific set of magnetotail parameters, the electron 
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finite Larmor radius and the electron drift term are respon- 
sible for the formation of a double peak in the total current 
density, even in those case where the proton current den- 
sity does not display any bifurcated structure and without 
turbulence. 

A number of experimental studies show that in many 
cases the electrons are the main current carriers [Kistler et 
al, 2005; Runov et al, 2006; Israelevich et al, 2008]. How- 
ever, some Geotail observations show that in some periods 
the ion current prevails, while, in other periods, the electron 
current prevails [Asano et al, 2003]. Keeping in mind that 
under different magnetotail conditions the main current car- 
riers can be either ions or electrons, here we try to assess 
the relative contribution to the cross tail current of differ- 
ent ion species, like protons and oxygen ions. The principal 
motivation to study + ions dynamics is the observation 
of rather energetic population of such ions, that shows a 
drift motion from dawn to dusk in the magnetic tail. The 
same observations [Wilber et al, 2004; Kistler et al., 2005; 
Cai et al, 2008] have shown that the + behavior is very 
different from that of the protons H + , analyzed till now, be- 
cause of the different Larmor radius (po+ = ^Ph+ : with the 
same initial energy), so that nongyrotropic effects are much 
more important for heavy ions. Out of the neutral sheet + 
ions and protons perform circular orbits around the mag- 
netic field lines. Nevertheless, since their Larmor radius is 
greater than the half thickness of the current sheet, the + 
ions can reenter in the neutral sheet, performing meandering 
orbits. An interesting situation arises where there is a thin 
current sheet and the plasma sheet is dominated by iono- 
spheric oxygen, for example during magnetospheric active 
periods: these are the most likely cases to find the plasma 
sheet dynamics dominated by non-adiabatic ions, that may 
also carry a large fraction of the current. In this paper we 
investigate the dynamics of ionospheric + ions in the near- 
Earth magnetotail in the presence and in the absence of 
magnetic turbulence. We find that + ions can support the 
formation of double peak even in the absence of magnetic 
fluctuations, as well as for relatively large values. 



the neutral sheet 11 times; at 09:48UT and at 09:50UT the 
current density showed a bifurcated structure [Runov et al., 
2003a] . At the same time a declining trend in proton density 
from 0.3 to 0.02 cm~ 3 is observed, that results in + domi- 
nance (n D + ~ 0.07cm -3 ) [Wygant et al, 2005; Wilber et al., 
2004]. This situation is rather common: Korth et al. [2002] 
have shown that + /H + energy density ratio of storm-time 
substorms exceeds 100%; for non-storm substorms the ratio 
is less and amounts to (15 — 65)%. Magnetic field value in 
the magnetospheric lobes and normal component of mag- 
netic field are approximately 20 nT and 5 nT. We have 
also estimated a magnetic fluctuations level <5B of around 
10 nT. The current sheet half-thickness measures approxi- 
mately 2500 km [Kistler et al, 2005]. 

The event of 17 August 2001 is another example in which 
the Cluster spacecraft were in the plasma sheet during a 
geomagnetic storm and observed the + outflow from the 
ionosphere into the tail [Sauvaud et al, 2004; Echer et al., 
2008]. In addition, during the storm main-phase, a sub- 
storm occurred. Oxygen beams were present in both the 
lobe and plasma sheet prior to this substorm, so this is an 
event where the plasma sheet is oxygen rich prior the sub- 
storm onset. After the substorm onset at 16:00 UT, the H + 
density and pressure decreased significantly, while + pres- 
sure remained high. From 16:36 to 16:46 UT, + was the 
dominant ion in the plasma sheet and showed a duskward 
motion. 

During the Cluster spacecraft crossing on 15 September 
2001, both an embedded proton (that manifests a pressure 
anisotropy mainly with py > p±) and a bifurcated oxygen 
ion thin current sheet (that exhibits a pressure anisotropy 
mainly with pu < p± and nongyrotropy) were observed 
[Cat et al, 2008]. Indeed, in the time interval (04:57:45- 
05:00:25)UT, a localized self-consistent current sheet equi- 
librium of oxygen ions was observed by CI and C4, althought 
in this very thin sheet the current contribution from oxygen 
ions is minor. Magnetic field value in the magnetospheric 
lobes and normal component of magnetic field are approxi- 
mately 25 nT and 3 nT. The estimation of the oxygen thin 
current sheet thickness is about 2500 km. 



2. Overview of oxygen observations 

Measurements have confirmed the presence of out-flowing 
ions + , H + , He + and other ions in every magnetospheric 
region; these ions have low energy in the ionosphere and 
high energy in the magnetosphere, thanks to several accel- 
eration mechanisms. It is known that the outflow of heavy 
ions from the ionosphere is strongly dependent on geomag- 
netic and substorms activity [Yau et al, 1985; Wilson et al., 
2004]. In particular + ions are not just significant in the 
inner magnetosphere, but actually the dominant ion species 
during magnetic storms [Daglis, 1991; Kistler et al., 2005]: 
they are observed to stream from dawn to dusk across the 
tail, carrying about 5 — 10% of the cross-tail current. Indeed, 
during storm times, + can dominate both the pressure and 
the density in the plasma sheet, that is already oxygen-rich 
because of the contribution from ion-outflow. 

Inspection of Cluster data shows that often the formation 
of the double current sheet is in association to the increase 
of the oxygen amount in the near magnetotail, for example 
during the events on 1 October 2001 [Runov et al., 2003a; 
Wilber et al, 2004; Wygant et al, 2005; Kistler et al, 2005], 
17 August 2001 [Sauvaud et al, 2004; Echer et al, 2008] 
and 15 September 2001 [Cai et al, 2008]. Generally, these 
layers are very thin, with typical thickness ranging between 
2000 and 5000 km. 

On 1 October 2001 at (0947-0951)UT, during a storm- 
time substorm, Cluster was located at [—16.5; 8.0; 0.5] Re 
(GSM). During this interval the Cluster barycenter crossed 



3. Test particle simulation 

The magnetic field model used in the simulation is de- 
scribed in full detail in Veltri et al. [1998]; Greco et al. 
[2002]. We recall here some of the main features. The 
considered magnetic field model consists of three terms, 
B = Boi(z) + B n + <5B(r). Ho x (z) is an unperturbed, sign 
reversing component, directed along the Earth-Sun axis. We 
use the expression of a modified Harris magnetic field rever- 

Bq = B tanh (f) -(f) cosh' 2 (^) 
°tanh(A)_(^)cosh- 2 (A) 

where Bo is the constant magnetic field in the magneto- 
spheric lobes, A = 0.251/ is the current sheet half thick- 
ness and L is the thickness of the considered magnetic field 
configuration (i.e. of the simulation box). B n = B n e z is 
the unperturbed, constant, normal component. <5B(r) cor- 
responds to magnetic fluctuations, here represented as the 
sum of static magnetic perturbations: 

SB{r) = <5B(k)e ff (k) exp [i (k • r + <&)], 

k,a 

where e a (k) are the polarization unit vectors, 0£ are ran- 
dom phases and 
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Table 1. Normalization quantities. 

Physical quantity Normalization Typic value 



Length 
Time 

Electric field 
Magnetic field 
Current density 
Velocity 
Temperature 
Ion density 



E 
Bo 
B /noL 
r E = cE /B 
m v V^/k B 



4000 km 

2 s 

0.2 mV/m 
20 nT 
4 nA/m 2 
10 km/s 
1.2 X 10 4 K 



cBq 1 4-KeLVE 0.5 cm 



where C is a normalization constant and l x , l y and l z are the 
turbulence correlation lengths in the x, y and z directions, 
respectively. The correlation lengths are fixed in terms of 
the thickness of simulation box as l x = l v = 0.25L and 
l z — 0.05L, in order to mimic the geometry of the mag- 
netotail (i.e. l z 3> l x ,l y ). The spectral index a for the 
near-Earth magnetotail is chosen as a = 2.3 [e.g. Zimbardo 
et al, 2009]. 

In addition to the above magnetic field configuration, we 
consider a constant cross tail electric field in the dawn to 
dusk direction, E y = E y e y , that cannot be removed by 
transforming to the de Hoffman- Teller frame for the pres- 
ence of three-dimensional magnetic fluctuations. Since we 
are considering static magnetic perturbations, the fluctuat- 
ing electric field is not included in the present runs. We may 
show that the fluctuating eletric field has a minor inpact on 
particle dynamics, because the ratio Va/vhi <C 1 [Veltri et 
al, 1998; Greco et al, 2002]. A study of time dependent 
fluctuations and the associated energization has been car- 
ried out by Greco et al. [2009]; Pern et al. [2009]. Instead 
in this study we are interested in studying the particles dy- 
namic and its role in the formation of the double current 
layer, rather than the acceleration mechanism to which these 
particles could be subject. In Table 1, normalizations and 
typical values of the physical quantities, inferred from ob- 
servations [Runov et al, 2006; Wygant et al., 2005; Kistler 
et al., 2005], are shown. 

We also denned SB = ^{SB ■ SB), with the average made 
over the simulation box, and 5b — SB /Bo and b n — B n /Bo- 
We consider that the source of the particles that are en- 
tering the current sheet from the lobes, at z = ±0.5L, is 
located somewhere in the magnetospheric mantle. The rela- 
tively cold ion distribution in the mantle magnetic field can 
be described as a shifted Maxwellian (e.g., [A shour- A bdalla 
et al, 1994]): 



they can carry a large fraction of the current away from the 
center z — 0. The nonzero B n plays the role of a guid- 
ing channel for the particles and causes them to exit the 
quasineutral sheet in the z direction [Speiser, 1965]. Typi- 
cally the particles, after traversing variable distances in y, 
stop meandering and travel towards the lobes with larger 
Larmor radius, since they gained energy being accelerated 
by E H (this effect is more apparent in the panels B-D). 

Increasing the values of magnetic fluctuations 8b (panels 
(B), (C) and (D)), the performed orbits become are per- 
turbed and the trajectories become more intricate, espe- 
cially for H + ions (see the blow up in Figure 1). In the 
regions of meandering motion, different deflections are ob- 
served, which slow down the motion along y and make the 
trajectory more tangled, although this effect is more evi- 
dent for protons. Starting with the same initial conditions, 
the + ions interact less with the magnetic turbulence than 
the protons H + , because of their greater Larmor radius: for 
T = 0.5 keV and for B = O.lBo (e.g. in the neutral plane), 
p(H+) ~ 1130 km and p(0 + ) ~ 4600 km. The turbulence 
characteristic lengths are \ m ax = L/k m in = 4000 km and 
^min = L/kmax = 333 km, where k max = 12 and k m i n = 1 
are the maximum and minimum wave number in each di- 
rection, respectively. The + ions are not influenced by the 
magnetic turbulence so much as protons; besides oxygen 
ions are able to transport current out of the neutral sheet, 
because of the amplitude of performed orbits. A similar dif- 
ference in the proton and oxygen interaction with magnetic 
turbulence was reported by Taktakishvili et al. [2007], with 
regard to plasma transport across the magnetopause current 
sheet. 

In order to gain understanding in the ion dynamics, we 
inject many particles in the simulation box and numerically 
integrate their equation of motion. Then we compute the 
distribution function moments as density n, current density 
j and temperature T on a three-dimensional grid with 1 grid 
point in x (many grid points in x could be used; however 
for the present runs we set this to 1, since a; is a statistically 
ignorable coordinate in our model), 601 in y and 40 in z. 
Then the distribution function moments have been averaged 
over y, the dependence on which is usually weak, to show the 
characteristic dependence on the z coordinate. In a test par- 
ticle simulation the normalization for the number density is, 
to good extent, arbitrary. Here it is based on the considera- 
tion that the ion current I y has to be strong enough to repro- 
duce the unperturbed magnetic field Bq x (z) introduced in 



/(«y,«x) = (\/2^4) 



exp 



^ u th 



where u is the streaming velocity along the unperturbed 
magnetic field and v t h is the thermal velocity. 

In the present simulation we assume b n = 0.1 for the 
normal magnetic field. Typically, 500000 particles are in- 
jected for each run with temperature Ti = 0.5 keV, in agree- 
ment with the values observed in the magnetospheric lobes. 
The corresponding ion injection velocities are respectively, 
vth{0 + )= 56 km/s and vth(R + )— 224 km/s [Vaisberg et al, 
1996]. Besides, we assume u — v t h for both species. In 
Figure 2, we show the projection on the yz plane of four 
sample trajectories of H + and + ions, for different values 
of magnetic fluctuations 5b. Once entered in the simulation 
box, in total absence of magnetic fluctuations (panel (A)), 
both types of particles move toward the central region z = 0, 
following the magnetic field lines while being subjected to 
the (E x B) drift. Inside of the quasi-neutral sheet, parti- 
cles begin to perform the typical meandering orbits, under 
the action of the unperturbed magnetic field, while being 
accelerated by the electric field E H . It's clear that + has 
larger Larmor radius than protons of the same energy and 




Figure 1. Blow up of + ions (black line) and pro- 
tons (red line) trajectories for b n = 0.1, E y = lEo and 
56 = 0.2. 
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Figure 2. Projection on the plane of four sample trajectories of + ions (black line) and protons (red 
line) for b n = 0.1, E y = 1E and Sb = (panel (A)), 8b = 0.1 (panel (B)), 8b = 0.2 (panel (C)), Sb = 0.6 
(panel (D)). 



our model. Indeed, Ampere's law leads to 2BoL x — 4nl y /c, 
with I y — J jydxdz the total current across a section of the 
current sheet of length L x . On the other hand, velocity is 
expressed in units of Ve, so that the normalization for den- 
sity is obtained an n* = Iy/[ e J Vydxdz], with I y constant 
for all the runs (this implies that, when the average value of 
V y is large, the density is low). 

In figures 3 and 4, we report the vertical profiles of oxygen 
temperature T and current density J y , for different values 
of electric field E y = (1, 2, 3, 4)E [Cattell and Mozer, 1982], 
magnetic fluctuations Sb = 0.0 (black line), 0.1 (red line), 
0.2 (blue line), 0.4 (purple line), 0.6 (green line) and for 
b n — 0.1, wich corresponds to a variety of observations in 
the near-Earth magnetotail. In figures 6 and 7, we show the 
vertical profiles of proton temperature T and current density 
J y for E y — (1, 2, 3, 4)Eq, b n — 0.1 and only for two values 
of the magnetic fluctuations Sb — 0.0 (black line), 0.6 (green 
line) (results for more values of Sb have been given by Greco 
et al. [2002]). Given the abundance of previous work on pro- 
tons, we restricted to show for this specie only run in which 
the oxygen shows a bifurcated current sheet, to be able to 
get the total current. The oxygen temperature in figure 3 in- 
creases with the electric field, but there is no clear influence 
of magnetic fluctuations on heating. Indeed, in each panel, 
different colour lines are at same level. A different situa- 
tion arises with protons. In this case, turbulence succeeds 
in scrambling the ordered velocity gained because of E y into 
all directions, leading to an effective thermalization of the 
potential drop [Greco et al., 2002]. The temperature profiles 
are larger in the center of the simulation box, where the en- 
ergization from the electric field is favored and the magnetic 
fluctuations are stronger. Finally, oxygen and proton tem- 
peratures are of the same order (~ 3 — 5 keV) in the case 
of E y = 4£"o and Sb — 0.6 (green line), in agreement with 
observations. 

From figure 4, it is clear that oxygen current profile is 
peaked and thin for zero to low fluctuation levels, and be- 



comes progressively broader as the fluctuation level is in- 
creased, while, at the same time, the maximum current value 
decreases. An obvious effect of magnetic fluctuations is, in- 
deed, to reduce the bulk velocity (and current density), by 
inducing random motions, which spread all around single 
particle velocities. Besides, oxygen current density assumes 
values which grow with the electric field. In the absence of 
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Figure 3. Vertical profiles for + ion temperature for 
Ey = (l,2,3,4)£o and b n = 0.1. Here 8b = 0.0 (black 
line), Sb = 0.1 (red line) e Sb = 0.2 (blue line), Sb = 0.4 
(purple line), Sb — 0.6 (green line). Temperatures are in 
units of m p 1/J/fc_B. 
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Figure 4. Vertical profiles for + ions current density for E y = (1, 2, 3, 4)Eo and b n — 0.1. Here Sb — 0.0 
(black line), Sb = 0.1 (red line) e 8b = 0.2 (blue line), 8b = 0.4 (purple line), 8b = 0.6 (green line). 




Figure 5. Total current density profiles for E y = (1, 2, 3, 4)i5o, b n — O.I.Here Sb — 0.0 (black line) and 
Sb = 0.6 (green line). The right column is for k = 5, the left one is for k = 2 (see text). 
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due to the velocity increase, caused by the electric field ac- 
celeration, so that the oxygen gyroradius and the effect of 
meandering orbits increases. The double peak disappears 
for 5b = 0.05 (not shown), 0.1 (red line), 0.2 (blue line), 
and appears again for levels 5b > 0.4 (purple and green 
line). The proton behaviour is similar in some aspects and 
opposite in other ones. We can see from figure 7 that the 
proton current density is larger than that of the oxygen one 
and, in contrast with the preceding case, it decreases with 
increasing electric field. This is due to the presence of the 
negative wings, that are more evident for small values of 
the electric field: indeed, since the integral I y of the ion cur- 
rent must remain constant, if the contribution of negative 
wings is big, the positive current must increase. A small 
hint of a double peak, in absence of magnetic fluctuations 
(black line), is present only for high values of electric field; 
instead, in presence of high magnetic fluctuations (5b — 0.6 
(green line)), protons support the formation of double layer 
for all E y values, the bifurcation being more evident with 
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Figure 6. Vertical profiles for H + ion temperature for 
E y = (1,2,3,4)£ and b„ = 0.1. Here 5b = 0.0 (black 
line) and 5b = 0.6 (green line). Temperatures are in units 
of m p V£/k B - 
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increasing of the electric field. Moreover, negative values 
in the current density profiles are observed: they are the 
so called diamagnetic currents, due to the magnetization 
current — cV x (-P±y^y), proportional to v\ [Zelenyi et ai, 
2000]. If the anisotropy is weak (v« ~ v±), particles will 
give a substantial contribution to the formation of the dia- 
magnetic wings; conversely, if there's a strong anisotropy 
(v\\ 3> v±), this contribution is small [Greco et ai, 2002; 
Zimbardo et ai, 2004]. 

In order to assess the oxygen contribution in the forma- 
tion of a bifurcated current sheet in the total current den- 
sity, we define the density ratio k = n(0 + )/n(H + ). Then, 
we use two different k values, wich are close to the real val- 
ues observed in the near-Earth magnetotail, when the oxy- 
gen is the most abundant species, to find the total current 
density. If k = 2, J tot = (2/3)J D+ + (1/3)J P ; if k = 5, 
Jtot = (5/6)J + + (1/6)J P . As we can see in figure 5 , we 
found that the double peak is present in the total current 
density profiles, also when it is not present in the proton 
current density, that is for E y = 2Eo and 5b = (black 
line). 

4. Discussion and Conclusions 

Motivated by the many observations of ionospheric oxy- 
gen + in the Earth's magnetotail, in this paper we have 
considered the dynamics of protons and oxygen ions in the 
magnetotail current sheet, both in the presence and in the 
absence of magnetic turbulence. The ion motion has been 
studied by a test particle simulation, where particles are in- 
jected in a magnetic quasi-neutral sheet with superimposed 
magnetic fluctuations. Several runs have been made in or- 
der to understand the differences between the proton and 
the oxygen behavior, with emphasis of the possibility to ob- 
tain a bifurcated current sheet. 

When protons are injected in the simulation box, a clear 
bifurcation of the current profile is obtained for turbulence 
5 levels SB/B > 0.3 (see also Greco et at [2002]). Con- 
versely, because of the non adiabatic meandering motion, 
oxygen ions are able to support a bifurcated current sheet 
also in the absence of magnetic fluctuations. One of the 
most intriguing results concerning the oxygen behaviour is 
the appearance of a double hump in the current profile both 
for the case of no magnetic turbulence and for relatively 
high levels of magnetic fluctuations. The explanation for 
this behaviour can be found in the interaction between oxy- 
gen ions and the current sheet. Looking at figure 8, which 
shows a typical oxygen trajectory for 5b = and E y — lEo, 
we can notice that orbit is very smooth and that oxygen 
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Figure 7. Vertical profiles for H + ions current density 
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ion probes regions of the current sheet far from the neu- 
tral plane (because of its large Larmor radius) during its 
meandering motion along y. Therefore, this effect is essen- 
tially related to the non adiabatic motion of oxygen ions. If 
a great number of this kind of trajectories are statistically 
added, there will be a concentration of particles with high 
v y at those distances, obtaining a current density J y which 
displays two peaks away from the neutral sheet [Delcourt 
and Belmont, 1998]. The distance between the two peaks 
should be of the same order of the typical excursion in z, 
that is d ~ ^/poX, where po is the oxygen Larmor radius in 
lobe magnetic field Bo [Greco et al, 2002]. From figure 4, 
the separation between the two bumps for the case 8b — 
(black lines) is ~ 0.2L = 800 km. If we consider that a typ- 
ical oxygen Larmor radius, in a magnetic field of 20 nT and 
with a temperature of 0.5 keV, is of the order of 800 km, we 
obtain for d a value of the order of 900 km, which is compa- 
rable with the distance between the two humps computed in 
the vertical profile of the current density. The influence of 
magnetic turbulence on J y is twofold: on the one hand, rela- 
tively low levels of magnetic perturbations (8b < 0.4) cause 
the two peaks to decrease and smear because fluctuations 
scatter the oxygen ions around. On the other hand, when 
the level of magnetic perturbations is high enough, the pres- 
ence of fluctuations bends the field lines also in y direction, 
and allows cross field motion away from the central plane. 

Another important difference between the proton and the 
oxygen behavior is found in the temperature profiles. While 
the proton temperature grows with the magnetic turbulence 
level, the oxygen temperature does not. More precisely, the 
oxygen temperature appears to be uniformly large (that is, 
independently of SB /Bo), in the sense that most of the po- 
tential drop, due to the electric field, is transformed into 
heat. This means that thanks to their larger Larmor ra- 
dius, + ions are able to gain a large fraction of the poten- 
tial drop. This is converted into disordered motion by the 
chaotic dynamics in the quasi-neutral sheet, while the influ- 
ence of turbulence is evident mostly in the outer part of the 
simulation box, see 3. Once again, the comparison of proton 
and oxygen dynamics shows that the former interact more 
with turbulence than the latter. To some degree, we can 
say that the turbulent fluctuations are averaged along the 
gyroorbit, a behaviour wich is also found in the studies of 
plasma transport in the presence of turbulence (e.g., Zim- 
bardo et al. [2006]; Pommois et al. [2007]). Finally, when 
realistic ratios of oxygen to proton are assumed, it is found 
that a bifurcated current sheets is obtained even in absence 
of magnetic fluctuations. This suggests that the presence of 
oxygen ions as the main current carriers can help to explain 
the observations of bifurcated current sheets in the magne- 
totail. On the other hand, a complete description of the 
magnetotail current should take into account the electron 
contribution (e.g., Greco et al. [2007]), as we plan to do in a 
future work. 
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